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The effect of geometrical confinement, atomic position and orientation of Silicon 
nanowires (SiNWs) on their thermal properties are investigated using the phonon 
dispersion obtained using a Modified Valence Force Field (MVFF) model. The spe- 
cific heat (C v ) and the ballistic thermal conductance (ftf a ') shows anisotropic variation 
with changing cross-section shape and size of the SiNWs. The C v increases with de- 
creasing cross-section size for all the wires. The triangular wires show the largest C v 
due to their highest surface-to- volume ratio. The square wires with [110] orientation 
show the maximum k\ since they have the highest number of conducting phonon 
modes. At the nano-scale a universal scaling law for both C v and K^ al are obtained 
with respect to the number of atoms in the unit cell. This scaling is independent 
of the shape, size and orientation of the SiNWs revealing a direct correlation of the 
lattice thermal properties to the atomistic properties of the nanowires. Thus, en- 
gineering the SiNW cross-section shape, size and orientation open up new ways of 
tuning the thermal properties at the nanometer regime. 
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I. INTRODUCTION 



The increasing variety of application of Silicon nanowires (SiNWs) ranging from MOSFETs^ 
to rechargeable batteries^, thermoelectric (TE) devices^, bio-sensors^, solar cells^, etc., re- 
quires a good understanding of the lattice thermal properties. Thermal requirements can 
be very contrasting depending on the type of applications, for eg. MOSFETs will require 
heat dissipation for better performance whereas thermoelectric devices require low thermal 
conductivity to maintain a good temperature gradient for higher TE efficiency^. 

Improvement in the process technologies have led to the fabrication of SiNWs with differ- 
ent shapes, sizes and channel orientations^^. SiNWs with channel orientations along [100] , 
[110] and [111] are the most commonly manufactured ones. At the nanometer scale, strong 
geometrical confinement, atomic positions and increased surface-to- volume ratio (SVR) play 
significant roles in determining the thermal properties of the SiNWs^^. 

Thermal transport measurements using techniques like the 3oo method^ and thermo- 
reflectance^have led to a good understanding of the heat flow in large nanowires (diameter > 
30nm). However, the physics of heat flow in ultra-scaled SiNWs is still not well understood^. 
The phonon spectra of SiNWs can provide a theoretical estimate of heat transport through 
ultra-scaled structured since phonons (lattice vibrations) are responsible for carrying most 
of the heat in semiconductors^* 10 * 11 * 14 ! 

In this paper we theoretically explore the effect of (i) cross-section geometry, (ii) cross- 
section size, and (iii) wire orientation of ultra-scaled SiNWs on the thermal properties such 
as the ballistic thermal conductance (k\ ) and the specific heat (C v ). Furthermore, analyt- 
ical expressions for the variation of these physical quantities with size for each cross-section 
shape and orientation are provided to allow for a compact modeling representation of ther- 
mal properties to be used in simulators like Thermal-Spice^ and Themoelectric module 
simulator^. 

The thermal properties of SiNWs differ considerably from bulk Si 6 - . The transition from 
the particle to wave nature of heat transport with structural miniaturization calls for im- 
proved heat transport models. The traditional continuum models for thermal conductivity 
by Callaway^ and Holland^ are based on the Debye limit for phonons, sound velocity and 
many other parameters, which render these models quite cumbersome at the nanoscale, as 
discussed by Mingo et at^. The dependence of continuum models on a large set of fit- 

2 



ting parameters make them unsuitable for predicting the lattice thermal properties with 
dimensional scaling. This limitation can be overcome by atomistic models which automati- 
cally take into account the effects of structural miniaturization like geometrical confinement, 
orientation effects, cross-sectional shape and surface-to-volume ratio effects^ -^. 

Theoretical estimates of the thermal conductivity and the specific heat using phonon 
spectra have been addressed in the literature in the past. Calculations of thermal conduc- 
tivity using the full phonon spectra of large to medium sized SiNWs along certain specific 
orientations have been performe d 10 * 19 ! The influence of surface roughness on nanowires has 
been studied in RefsP^Zl. A study of specific heat in [111] SiNWs is reported in RefP^. 
Mingo et. al theoretically bench-marked the thermal conductivity of large SiNWs (diameter 
> 34nm)~° and also studied the effect of amorphous coating on the thermal conductivity 
of SiNWs 11 . A complete study of orientation effects on the thermal conductivity of SiNW 
has been provided in RefP. The thermal conductivity of Si nano-cluster d 24 * 25 -* and hollow Si 
nanowire d 14 * 26 ! has also been studied. 

In this work we utilize a semi-empirical phonon model to perform thermal calculations. 
The phonon dispersion is obtained using a Modified Valence Force Field (MVFF) model™ 1 . 
The continuum models^^l though computationally simpler and faster, lack the proper 
physics to extend them to nanostructures due to the use of open fitting parameters^. The 
first principles based models^QLiU are limited to extremely small structures (W<2nm) with 
few thousand atoms, and are computationally very expensive. The choice of the MVFF 
model is guided by the following factors, (i) an adequate amount of physics to understand 
the thermal properties, and (ii) a lower computational cost thus, allowing the model to 
handle up to few million atoms for phonon calculations in realistic structures^. At the 
same time the MVFF model also explains the experimental phonon and elastic properties 
in zinc-blende materials very wel l 27 1 33 ' ^ . 

The emphasis of this work is to purely understand the atomistic effects on the thermal 
properties of ultra-scaled SiNWs such as the effect of the surfaces and atomic coordination 
number. This understanding is further translated into analytical equations allowing easy 
usage for other thermal modeling like thermo-electric circuit simulations which require C v 
and Ki as inpulP^. Scattering effects are dominating at the nanometer scale in determining 
the lattice thermal conductivity. However, scattering effects have been neglected in the 
present study in order to understand the geometric and atomistic effects which will be 
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otherwise convoluted by scattering. Also the ballistic phonon transport can be supported 
in SiNWs with width smaller than 20nm as pointed in RefP^. The scattering effects could 
be studied as a future work. 

The paper is organized as follows. Section [TT] provides a brief description of the MVFF 
model for the calculation of phonons in SiNWs, the calculation of thermal properties 



(Sec. II B) and details of the SiNWs used for the study (Sec. II C). The effect of cross- 
section shape, size and wire orientation on the specific heat and thermal conductance are 



discussed in Sec. HI A and Sec. Ill B , respectively. This is followed by a discussion on the 



atomistic effects on the thermal properties in Sec. |III C| Conclusions are summarized in 
Sec. ED 



II. THEORY AND APPROACH 

A. MVFF Phonon model 

In the MVFF model, the phonon frequencies are calculated from the forces acting on 
atoms produced by finite displacements of the atoms from their equilibrium positions in a 
crystal. First the total potential energy of the solid (U) is estimated from the restoring 
force(F). In the MVFF model, U is approximated as^, 



r j^fc 

p»5_E E °£+ E MP 

iTfjik , jjjik \ , \ ^ jjjikl 
> u bs-bs u bs-bb) 2-^1 bb-bb 

j,k,l£COPi 



(1) 



where Na, nn(i), and COPi represent the total number of atoms in one unit cell, number 
of nearest neighbors for atom 'i', and the coplanar atom groups for atom 'i', respectively. 
The first two terms U l J s and Ul b tk represent the elastic energy obtained from bond stretching 
and bending between atoms connected to each othei^. The terms U bs _ hs , L^_ 66 , and U^_ bb 
represent the cross bond stretching^, the cross bond bending-stretchingP, and the coplanar 
bond bending^ interactions, respectively. The detailed procedure for obtaining the phonon 
spectra in bulk Si and NWs are outlined in RefP 33 . 
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B. Lattice thermal properties 

The complete phonon dispersion provides information about the thermal properties of 
nanostructures^^. The constant volume temperature (T) dependent specific heat (C V (T)) 
can be evaluated using the following relational, 



C V (T) = (2) 

kR r (^rL£) . e xp(^^) 1 
(-)TP^ kBT 1 [J/kgK], 

where ks, h and the Boltzmann's constant, reduced Planck's constant, and the 

mass of the SiNW unit cell in kg, respectively. The quantity u) n>q is the phonon frequency 
associated with the branch 'n' and crystal momentum vector 'q'. 

For a semiconductor slab/wire with a small temperature difference AT between its two 
extremities, the thermal conductance (fi^ aZ ) is obtained using Landauer's methocP^ a ij 6 | 10 | 1 H 



n b r(T) = h 



d_ 
df 



(exp( 



T(u) ■ M(uo) ■ u 



(3) 



k B T' 



1 



dco [W/K], 



The term M(u) is the number of phonon modes at frequency u and T(u) is the trans- 
mission for each mode. For ballistic conductance each mode transmits with a probability of 
1 while for coherent scattering dominated conductance the transmission value is less than 1. 



C. Si Nanowire details 



In this study, four types of cross-section shapes for [100] SiNWs have been considered 
namely, (a) circular, (b) hexagonal, (c) square and (d) triangular (Fig. [T]). Square SiNWs 
with [110] and [111] channel orientations have been studied too (Fig. [2]). The feature size 
is determined by the width parameter W. The value of W is varied from 2 to 6 nm. The 
confinement direction of the SiNWs are along Y and Z. The heat transport direction is 
along X. The surface atoms are allowed to vibrate freely without any passivating species. 
These wires are assumed to have a tetrahedral geometry. It has been shown that wires with 
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diameter below 2nm tend to lose the tetrahedral structur e 1 40 * 41 ! due to surface pressure and 
internal strain. 



III. RESULTS AND DISCUSSION 

In this section the results on the effect of cross-section shape, size and orientation on 
the thermal properties of SiNWs are presented and discussed. All the thermal quantities 
are calculated at 300K. However, the analysis holds for any temperature (T) where the 
anharmonic phonon effects are small. For T < Toebye the anharmonic interactions are quite 
smalP^. For bulk Si, Toebye ~725Kp2] anc l it further increases for SiNWs 321 . 

A. Specific heat in SiNWs 

Influence of the shape and size on C v : The C v of SiNWs increases with decreasing cross- 
section size in all the wire shapes (Fig. |3^i) . The size dependence of C v can be approximated 
by the following relational, 

C v (W) = C% M + £, (4) 

where, A is a fitting parameter extracted from the linear fit of the numerical simulations. 
Table [i] shows the value of C^ ulk and A for each geometry. As W — > oo (increasing cross- 
section size), the C v of all the SiNWs converges to a fixed value of ~681 J/kg.K which 
is reasonably close to the experimental C v value for bulk Si (~682 J/Kg.K as provided in 

Refs [23]l3]) triangular 

wires show the maximum C v for all the W value, whereas the 
other shapes show similar C v values at any given cross-section size (Fig. |3^i). 

The plot of AC V {=Cf Te - C^ lk ) vs. SVR (SVR = Total surface Atoms/Total atoms in 
unit cell) shows a linear behavior (Fig. [3]d), which can be represented as, 

C7 re « m c x SVR + C b v ulk , (5) 

where m c describes the additional contribution to the C v of the SiNWs with increasing 
surface-to-volume ratio. The value of m c is positive for all the wire shapes (Fig. [3|d) which 
corroborates the fact that specific heat increases with increasing surface areaPl Different 
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coordination number of surface atoms for the various cross-section shapes result in different 
m c values which depict the atomistic effect on the C v value in ultra-scaled SiNWs. 

The C v increase with decreasing W can be attributed to two phenomena, (i) phonon 
confinement due to small cross-section size and (ii) an increased surface-to-volume ratio 



energy bands active at a given temperature (see Eq. |2|. Thus, more energy is needed to 
raise the temperature of the smaller wires. 

The shape dependence of the C v can be understood from Eq. [5j The variation of (i) SVR, 
and (ii) m c with W for different shapes decide the eventual C v order. Figure |4]^a) shows 
that triangular wires have the maximum SVR while the other shapes have similar SVR at a 
fixed W. The increasing SVR results in a higher phonon density of states (DOS) associated 
with the wire surface, which further enhances the specific heat with decreasing wire cross- 
section 23 . The square wires provide the largest surface contribution to C v as depicted by 
the variation in m c (AC^/SVR) (Fig. |4]d). An optimal value of SVR and m c in SiNWs 
will maximize the C v . The SVR x m c value has the following order, triangle (~36) > 
hexagonal (~29) > square (~27) > circular (~26). Thus, triangular wires depict the highest 
C v due to the highest SVR. However, the trends for SVR (SVR sg < SVRh ex < SVR ci ) 
and m c (m s c q > m h c ex > m™) are opposite for the other shapes, hence resulting in almost 
similar C v values. The C v values have the following shape order in [100] SiNWs: triangular 
> hexagonal »s square »s circular. 

Influence of orientation on C v : The specific heat is also a function of the SiNW orientation 
(Fig. [5^i). The C v varies inversely with W, similar to the trend extracted from different 
shapes (Eq. [1]). The width parameters (Eq. [4]) for the variation of C v with orientation are 
provided in Table [TTJ The AC,, value again shows a linear variation with SVR for different 
wire orientations (Fig. [5Jd). The C v has the following trend with orientation for different 
cross-section sizes, C* 00 > C^ w > C* 11 . This trend can be explained again by looking 
at the impact of (i) SVR and (ii) m c ( Eq. [5]) on the overall C v value. The SVR shows the 
following order with W, SVR 111 > SVR 110 > SVR 100 (inc. of ~1.2x from [100] to [111]) as 
illustrated in Fig. [6^1. However, m c shows the following order with W, m^ 00 > mj 10 > mj 11 
(dec. of ~1.7x from [100] to [111]) (Fig. ^jp). The larger surface to volume ratio plays the 
main role in deciding the C v trend for SiNWs with different orientations. 
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B. Ballistic thermal conductance of SiNWs 



The thermal conductance indicates how a structure can carry heat. For high thermo- 
electric efficiency (ZT) a small thermal conductance is needed ^ 11 * 44 * whereas for CMOS 
transistors a high k\ is preferred to evacuate the heat. The cross-section shape, size and 
wire orientation of SiNWs can be used to tune their the thermal conductance. 

Influence of the shape and size on n\ al : The ballistic thermal conductance (/Cj ) for 4 
different shapes (Fig. [TJ of [100] SiNWs are calculated as shown in Fig. [7^. The K^ al can be 
fitted according to the following size (W) relation, 



, , ,'W\ d 
K b l al (W) = K (-) , (6) 
ao / 



where ao is the silicon lattice constant (0.5431 nm), d is a power exponent, and k,q is a 
constant of proportionality. The values of 'd' and k for different wires shapes are provided 



in Table III The value of d varies between 1.92 and 2.011 which implies that «j has a 
similar size dependence for all the wire shapes. However, the pre-factor value (kq) reflects 
the shape dependence. This pre-factor has the same ordering as the thermal conductance 



ordering (Fig. [7^ and Table III). 



The ballisitic thermal conductance exhibits a linear behavior with the number of atoms 
per unit cell (NA) as depicted in Fig. This linear relation can be approximated by the 
following equation, 

k^ 1 ^m k x NA + ^{0), (7) 

where the slope represents the average contribution from each atom in the unit cell to 
Ki al and Kj a '(0) is the thermal conductance at NA = 0. The value of «j (0) is zero within 
numerical error (K bal (0) ~ le — 7) which is expected for NA = 0. The value of m k takes into 
account the surface, shape and atomic effects since the calculation procedure involves the 
complete phonon dispersion. This relation shows a direct correlation of the atomistic effects 
to the ballistic thermal conductance. 

The size dependence can be explained by the fact that the larger wires have (i) more 
phonon sub-bands resulting in higher number of modes (M{uS) in Eq. [3]) and (ii) a higher 
acoustic sound velocity which is responsible for a larger heat conduction in these SiNWs^. 

The shape dependence can be explained as an interplay of two effects, (i) the total number 
of atoms (NA) present in the unit cell of SiNWs, and (ii) the average contribution of every 
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atom towards k\ . For a fixed cross-section size W, the NA ordering is NA sq > NA& > 
NAhex > NATeri (Fig- [8^)- The total number of phonon branches are 3xNA, due to 
the three degrees of freedom associated with each atompSl. The number of phonon modes 
(M(u)) is directly proportional to the phonon energy sub-bands. The contribution per atom 
to the thermal conductance (rrik) stays almost constant with the wire cross-section size for 
a given shape (Fig. [8b). Since, the values of are quite similar for all the cross-section 



shapes , the ordering of NA with shape governs the dependence of K b on the cross-section 
shape. 

Influence of wire orientation on K bal : The thermal conductance is anisotropic in SiNWs 
with the following order, k] w > nj 00 > nj 11 (Fig. [9^). This result is similar to the one 
reported in RefP. The Ki value exhibits a linear variation with NA for all the wire orientations 
(Fig. |9}3). The width parameters for the thermal conductance (Eq. [6] ) for different wire 



orientations are provided in Table IV The order of the ballistic thermal conductance with 



W for different orientations can be understood by the product of NAxm& (P nm ). The NA 



has the following variation, NA 111 > NA 110 NA 100 (Fig. 10 1), whereas depicts the 



following order m\ 00 > > ml 10 > ml 11 (Fig, 10 d). These two orders are opposite to each 
other. However, the product shows the following order, > Pnm > Pnm)- Thus, [110] 

wires give the highest K\ due to the optimal value of NA and m^. 

An important point to note is that K b t al is expected to decrease further in smaller wires due 
to phonon scattering by other phonons, interfaces and boundarie j 10 * 11 ^ which are neglected in 
this present study. The main idea here is to understand the geometrical effects on the phonon 
dispersion and the lattice thermal properties of these small nanowires which is attributed 
to (i) the modification of the phonon dispersion, and (ii) phonon confinement effects in the 
coherent phonon transport regime. 

C. Discussion 

In this work, all the thermal properties are shown to scale with W. Also NA depends on 
W as follows, 



NA oc W\ 

where 7 > 0. So using Eq. (§, Q and @ can be recasted in terms of NA as, 
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(8) 



AC V (NA) = C -(NA)^ (9) 

= C • (NA)-** 
nf d (NA) = K ■ (NA)y (10) 
= K ■ (NA)", 

where, Cq and K are the pre-factors. Thus, a universal power law can be derived for the 
thermal properties depending on the number of atoms per unit cell (NA) which represents 
the atomistic effect on the thermal quantities. In these SiNWs, 1.98 < 7 < 2.1 which 
provides the limits for p and 77, 



p E [0.93 , 1.03] (11) 
77 E [0.48 , 0.50] (12) 

The variation in the thermal conductance with NA is illustrated on a log-log scale in Fig. 



111. All the SiNWs depict almost the same power law with an average exponent value of 



0.97, which is in the limit derived in Eq. 11 Similarly the variation in C v with NA is 



plotted on a log-log scale in Fig. [TTJd. All the SiNWs used in this study show the same 
power law (Eq. [9]) with an average exponent of -0.51, which is in the limit derived in Eq. 



12 Thus, the thermal quantities show a universal power law behavior with the number of 



atoms in the unit cell (NA) irrespective of the details of the unit cell. The details of shape 



and orientation are embedded in the pre-factors Co and Kq (Eq. [9] and 10). 

The relationship of C v and Ki to the shape, size and orientation of SiNWs have been 
provided explicitly in Eqs. Q, ([5]), (J6J), and ([7]). These closed form analytical expressions 
are very handy for the compact modeling of the thermal and thermoelectric properties of 
SiNWs™ 1 . S ince these expression are derived from physics-based model, they capture the 
important geometrical and atomistic effects, thus enabling fast modeling of realistic systems. 

IV. CONCLUSIONS 



We have shown the application of the MVFF model for the calculation of the thermal 
properties of SiNWs. It has been shown that at the nanometer scale these thermal prop- 
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erties are quite sensitive to the wire cross-section size, shape, and orientation. Analytical 
expressions for the size dependence of the thermal properties of SiNWs of different cross- 
section shape and channel orientation have been provided. They can be used as component 
for the compact modeling of the thermal properties of ultra-scaled SiNWs. It has been 
demonstrated that all the SiNWs follow a universal power law for the specific heat and the 
thermal conductance which reveals the impact of the atomistic details on these properties. 
The triangular SiNWs show a high C v and low Ki, thus making them good candidates for 
thermoelectric devices. The [110] oriented square Si nanowires are better in terms of heat 
dissipation due to their high thermal conductance and are therefore good candidates for 
transistors from a heat management point of view. 
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TABLES 



TABLE I. Width Parameter for C v in [100] SiNWs at T = 300K. 



Shape 


Qbulk 


A 


units — > 


(J /kg -K)(J- 


nm/kg • K) 


Circular 


681.4 


47.82 


Hexagon 


681.2 


53.58 


Square 


680.9 


52.66 


Triangular 


681.0 


73.43 



TABLE II. Width Parameter for C v in SiNWs at T = 300K. 



Orientation 




A 


units — > 


(J /kg -K) {J- 


nm/kg • K) 


[100] 


680.9 


52.66 


[110] 


681.5 


43.03 


[111] 


681.2 


39.4 



TABLE III. Width Parameter for n\ al in [100] SiNWs at T = 300K. 



Shape Ko 


(nW/K) 


d 


Circular 


0.133 


1.92 


Square 


0.141 


1.96 


Triangular 


0.062 


1.97 


Hexagon 


0.097 


2.01 



TABLE IV. Width Parameter for n\ al in SiNWs at T = 300K. 



Orientation kq 


(nW/K) 


d 


[100] 


0.141 


1.96 


[110] 


0.204 


1.87 


[111] 


0.129 


1.88 
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FIG. 1. Projected unit cell structures of free-standing [100] oriented silicon nanowires with (a) 
Circular, (b) Hexagonal, (c) Square, and (d) Triangular cross-section shapes. Width and height of 
the cross-section are defined using a single width variable W (width = height). These structures 
are at W = 2nm. 



[100] [110] [111] 




[a] [b] [c] 

FIG. 2. Projected unit cells of square free-standing SiNWs with (a) [100], (b) [110], and (c) [111] 
wire axis orientation. The width and the height of the SiNWs are defined using W. Here, W = 
2nm. 
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FIG. 3. (a) Dependence of the specific heat (C v ) on cross-section shape and size in [100] SiNWs. 
(b) Variation in AC V ( = C v - C b v ulk ) with SVR for all the [100] SiNW shapes. C b v ulk for each 
shape is given in Table. [Ij 
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FIG. 4. (a) Surface-to-volume ratio (SVR) for different cross-section shape and size [100] SiNWs. 
(b) Incremental contribution to the specific heat with SVR for different cross-section shape [100] 
SiNWs. 




FIG. 5. (a) Dependence of the specific heat (C v ) on the orientation of square SiNWs. (b) Variation 
in AC V ( = C v - C b v ulk ) with SVR for all the SiNW orientations. C b v ulk are taken from Table, g 
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FIG. 6. (a) Surface-to-volume ratio (SVR) for different wire orientations of square SiNWs. (b) 
Incremental contribution to the specific heat with SVR for different orientations of SiNWs. 
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FIG. 7. (a) Effect of cross-section shape and size on the ballistic thermal conductance of [100] 
SiNWs. (b) Variation in K^ al with the total number of atoms per unit cell (NA) for different 
cross-section shapes. 
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FIG. 8. (a) The number of atoms (NA) in [100] SiNW unit cell for different cross-section size and 
shapes, (b) Contribution to ki per atom for different cross-section shapes. 
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FIG. 9. (a) Effect of size on the ballistic thermal conductance in SiNWs with different orientations, 
(b) Variation in with the total number of atoms in the unit cell (NA) for different wire 
orientations. 
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FIG. 10. (a) The number of atoms (NA) with W in one SiNW unit cell for different orientations, 
(b) Contribution to ki per atom for different SiNW orientations. 
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FIG. 11. (a) The number of atoms (NA) in one [100] SiNW unit cell for different cross-section size 
and shapes, (b) Contribution to m per atom for different cross-section shapes. 
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